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Abstract 

I. J. Schoenberg proved that a function is positive definite in the unit 
sphere if and only if this function is a nonnegative linear combination 
of Gegenbauer polynomials. This fact play a crucial role in Delsarte's 
method for finding bounds for the density of sphere packings on spheres 
and Euclidean spaces. 

One of the most excited applications of Delsarte's method is a solution 
of the kissing number problem in dimensions 8 and 24. However, 8 and 
24 are the only dimensions in which this method gives a precise result. 
For other dimensions (for instance, three and four) the upper bounds 
exceed the lower. We have found an extension of the Delsarte method 
that allows to solve the kissing number problem (as well as the one-sided 
kissing number problem) in dimensions three and four. 

In this paper we also will discuss the maximal cardinalities of spherical 
two-distance sets. Using the so-called polynomial method and Delsarte's 
method these cardinalities can be determined for all dimensions n < 40. 

Recently, were found extensions of Schoenberg's theorem for multivari- 
ate positive-definite functions. Using these extensions and semidefinite 
programming can be improved some upper bounds for spherical codes. 



1 Introduction 

Let M be a metric space with a distance function r. A real continuous function 
f(t) is said to be positive definite (p.d.) in M if for arbitrary points p%, .. . ,p r 
in M, real variables x\, . . . , x r , and arbitrary r we have 

r 

<>3=1 

or equivalently, the matrix (/(ty)) ^ 0, where the sign >z stands for: "is 
positive semidefinite" . 

Let §™ _1 denote the unit sphere in K n , and let tf>ij denote the angular dis- 
tance between points Pi,Pj- Schoenberg [33] proved that: 

/(cos if) is p.d. in S™ _1 if and only if f(t) = X)?*=o fkG^ (t) with all ff. > 0. 
Here G^\t) are the Gegenbauer polynomials. 
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Schocnbcrg's theorem has been generalized by Bochner [H] to more general 
spaces. Namely, the following fact holds: / is p.d. in a 2-point-homogenous 
space M if and only if f(t) is a nonnegative linear combination of the zonal 
spherical functions $k(t) (see details in [TH1 Th. 2], [TU1 Chapter 9]). 

Note that the Bochner - Schoenberg theorem is widely used in coding theory 
and discrete geometry for finding bounds for error-correcting codes, constant 
weight codes, spherical codes, sphere packings and other packing problems in 
2-point-homogeneous spaces (see [TU1 HH [211 1251 123 132] and many others). 

The paper is organized as follows: 

Section 2 recalls definitions of Gegenbauer polynomials and considers Del- 
sarte's method for spherical codes. 

Section 3 discusses applications of Delsarte's method for the kissing number 
problem. One of the most excited applications of Delsarte's method is a solution 
of the kissing number problem in dimensions 8 and 24. However, 8 and 24 are 
the only dimensions in which this method gives a precise result. For other 
dimensions (for instance, three and four) the upper bounds exceed the lower. 
We have found an extension of the Delsarte method that allows to solve the 
kissing number problem (as well as the one-sided kissing number problem) in 
dimensions three and four. 

Section 4 discusses maximal cardinalities of spherical two-distance sets. Us- 
ing the so-called polynomial method and Delsarte's method these cardinalities 
can be determined for all dimensions n < 40. 

Section 5 considers Sylvester's theorem and semidefinite programming (SDP) 
bounds for codes. Delsarte's method and its extensions allow to consider the 
upper bound problem for codes in 2-point-homogeneous spaces as a linear pro- 
gramming problem with perhaps infinitely many variables, which are the dis- 
tance distribution. We show that using as variables power sums of distances 
this problem can be considered as a finite semidefinite programming problem. 
This method allows to improve some linear programming upper bounds. 

Section 6 discusses an application of the extended Schoenberg's theorem 
for multivariate Gegenbauer polynomials. This extension derives new positive 
semidefinite constraints for the distance distribution which can be applied for 
spherical codes. 

2 Delsarte's method 

2-A. The Gegenbauer polynomials. Let us recall definitions of Gegenbauer 
polynomials. Let polynomials C[ (t) are defined by the expansion 

(1 - 2rt + r 2 )< 2 -")/ 2 = jr r*c£ B) (t). 

fc=0 

Then the polynomials := C[ (t)/Ci (1) are called Gegenbauer or ul- 

traspherical polynomials. (So the normalization of is determined by the 
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condition Gj, (1) = 1.) Also the Gegenbauer polynomials G^f 1 can be defined 
by the recurrence formula: 

Mn)_i r {n)_, r (n)_ (2k + n - A) t G^, - (k - 1) C{% 

G -1, G x -i, ...,G fc - fc + „- 3 

Note that for any even fe > (resp. odd) G^ (t) is even (resp. odd). 
Therefore, G 2 k and G 2 f+i are orthogonal on [—1,1]. Moreover, all polynomials 
G^ are orthogonal on [—1,1] with respect to the weight function (1 — £ 2 )(™~ 3 )/ 2 : 

G { k n) (t) Gf ] (t) (1 - i 2 )(«- 3 )/ 2 dt = 0, ^ £ 

Recall the addition theorem for Gegenbauer polynomials: 
G^ ( cos $1 cos 02 + sin 0i sin 02 cos ip) 

k 

= $> nfcs Gtt 2s) (cos0i) G< n _+ 2s) (cos0 2 ) (sin^r (sin0 2 ) s G^ (cos p), 

s=0 

where c n /c S are positive coefficients whose values of no concern here (see [9l 114] V 

2-B. Schoenberg's theorem. Using the addition theorem for Gegenbauer 
polynomials Schoenberg |34j proved the following theorem: 

Theorem 2.1. A continuous real function f (cos (p) is positive definite in § ra_1 if 
and only if f is a nonnegative linear combination of the Gegenbauer polynomials, 
i.e. f{t) = £~ „ fkG { k n) (t) with all f k > 0. 

This theorem also can be proved using the addition theorem for harmonic 
polynomials (see details in [32]). 

2-C. Delsarte's inequality. Let {p±,p2, ■ ■ ■ ,Pm} be any finite subset of the 
unit sphere §™ _1 . By ipij = dist(pi,pj) we denote the spherical (angular) dis- 
tance between pi, pj. Clearly, cos ipij = (pi,p 3 ). 

If a symmetric matrix is positive semidefinite, then the sum of all its entries 
is nonnegative. Schoenberg's theorem implies that the matrix (G^ (%)) is 
positive semidefinite, where tij :— cos <pij, Then 

M M 

EE G i n) fe)>0 (2-1) 

i=l j=l 

Suppose a continuous functions / : [— 1, 1] — * R is p.d. in S™ _1 . Then 

oo 

f(t) = J2hG[ n) (t) 

k=Q 
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with all fk > 0. Let 

M M 

S(X) = S f (X) : X» :/ <- : - 

i=l j=l 

Using (2.1), we get 

oo / M M \ M M 

fc=0 \i=l i=l / i=l i=l 

Thus 

S f (X) > f M 2 . (2.2) 

2-D. Delsarte's bound. We say that X = {pi, . . . ,pm} C S™^ 1 is a spherical 
ip-code, where < ip < w, if for all i 7^ j, Uj = cos cfoj < z := cos^, i.e. 
tjj 6 [— 1, z]. In other words, the angular separation between distinct points 
from X is at least ip. Denote by A(n, ip) the maximal size of a f/'-code in 

Theorem 2.2 ( [TTJ fTJl [TB] ) . Let a continuous function f : [—1, 1] — > R fee p.d. 
m S"" 1 . Let /(£) < /or at/ £ £ [-1, cob ip]. Then 

A(n,V)<4^- 
Jo 

Proof. Let X = {pi, . . . ,pm} C § n_1 be a spherical i/'-code. Clearly, /(£«) = 
/(l). By assumptions we have f(Uj) < for all i ^ j. Therefore 

S f (X) = M/(l) + 2/(ti 2 ) + . . . + 2/(£ M -i,m) < M/(l). 

If we combine this with (2.2), then we get M < /(l)//o- d 



3 The kissing problem 

3-A. The kissing number problem. The kissing number problem asks for 
the maximal number k(n) of equal size nonover lapping spheres in n-dimensional 
space that can touch another sphere of the same size. In other words, k{n) = 
A(n, 7r/3), i.e. k(n) is the maximal size of a spherical 7r/3-code of length (di- 
mension) n. 

This problem in dimension three was the subject of a famous discussion 
between Isaac Newton and David Gregory in 1694. In three dimensions the 
problem was finally solved only in 1953 by Schutte and van der Waerden [56] , 

In 1979 Levenshtein [TI5], and independently Odlyzko and Sloane [31] (= 
[TUl Chap. 13]), using Delsarte's method, have proved that fc(8) = 240, and 
fc(24) = 196560. Moreover, Bannai and Sloane [S] (= (TUl Chapter 14]) proved 
that the maximal kissing arrangements in these dimensions are unique up to 
isometry. However, n = 8, 24 are the only dimensions in which this method 
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k=Q 



gives a precise result. For other dimensions (for instance, n = 3, 4) the upper 
bounds exceed the lower. 

3-B. The kissing problem in dimensions 8 and 24. The proofs in [THl 
|3"T] that fc(8) = 240 and fc(24) = 196560 are surprisingly short, clean, and 
technically easier than all known proofs in three dimensions. Indeed, let 

Mt) = (t - 1/2) t\t + i/2) 2 ( t +i) = E /i 8) 4 8) (*)> 

and 

10 

Mt) = (f - l/2)(t 1/4) 2 t\t + l/4) 2 (i + l/2) 2 ( £ + 1) = £ fGf 4 ' (i). 

fe=0 

Since all /<; 8) > 0, /f 4) > and / 8 (i) < 0, / 24 (i) < for all t e [-1,1/2] 
Theorem 2.2 yields 

*(8) < - 240, fc(24) < - 196560. 

/o /o 

For rt = 8, 24 the minimal vectors in sphere packings E$ and Leech lattices give 
these kissing numbers. Thus fc(8) = 240, and fc(24) = 196560. 

3-C. The kissing problem in four dimensions. It is not hard to see that 
fc(4) > 24. Indeed, the unit sphere in R 4 centered at (0,0,0,0) has 24 unit 
spheres around it, centered at the points (±\/2> ±v2> 0, 0), with any choice of 
signs and any ordering of the coordinates. The convex hull of these 24 points 
yields a famous 4-dimensional regular polytope - the "24-cell" . Its facets are 24 
regular octahedra. 

Let f OS (t) = h+hGf\t)+. . .+f 9 G ( *\t), where / = 1, /i = 3.6181, / 2 = 
6.1156, h = 7.0393, / 4 = 5.0199, f 5 = 2.313, / 6 = f 7 = f 8 = 0, f 9 = 0.4525. 
This polynomial was applied by Odlyzko and Sloane [31] to prove that fc(4) < 25. 
Since fos(t) < for t £ [—1, 1/2], Delsarte's bound gives 

fc(4) = A(4,tt/3) < fos(l)/fo = fos(l) « 25.5584. 

Thus, fc(4) < 25. 

Note that Arestov and Babenko [I] proved that the bound fc(4) < 25 cannot 
be improved using Delsarte's method. 
Let 

1344 n 2688 7 1764 - 2048 4 1229 , 516 , 217 2 

fJt) := t 9 t 7 H t 5 H 1 4 t 3 t 2 t 

1 v ' 25 25 25 125 125 125 500 125 

In [26] we proved that fc(4) = 24. This proof is based on the following two 
lemmas: 
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Lemma 3.1. Let X = {x\, . . . , xm} be points in the unit sphere S 3 . Then 

M M 
i=l i=i 

Proof. The expansion of /4 in terms of [/& = is 

A- E />, + 2 * + f u 2 + |i V, + f * + 1 w . 

We see that all /f 4) > and / (4) = 1. So Lemma 3.1 follows from (2.2). □ 



Lemma 3.2. Suppose X — {x±, . . . , xm} is a subset o/S 3 such that the angular 
separation between any two distinct points at least tt/3. Then 

M M 
z=l 3=1 

It's not easy to prove this lemma. A proof is given in [2BJ Sections 4,5,6]. 
Theorem 3.1. fc(4) = 24. 

Proof. Let X be a spherical 7r/3-code in S 3 with M — fc(4) points. Then X 
satisfies the assumptions in Lemmas 3.1, 3.2. Therefore, M 2 < S(X) < 25M. 
From this M < 25 follows, i.e. M < 24. From the other side we have fc(4) > 24, 
showing that M = fc(4) = 24. □ 

3-D. The kissing problem in three dimensions. Our extension of the 
Delsarte method can be applied to other dimensions and spherical ip-codes. The 
most interesting application is a new proof for the Newton-Gregory problem, 
A:(3) < 13. In dimension three all computations are technically much more easier 
than for n = 4 (see [2"4]). 
Let 

, . , 2431 n 1287 7 18333 » 343 4 83 , 213 , t 1 

fJt) = 1 9 1 7 + 1 5 + 1 4 1 3 1 2 + . 

J v ' 80 20 400 40 10 100 10 200 

Then for any M-point kissing arrangement X we have Sf 3 (X) < 12.88M (see 
details in [531 HI])- The expansion of f% in terms of Legendre polynomials 
Pk = G ( k 3) is 

f 3 = P + I.6P1 + 3.48P 2 + 1.65P 3 + 1.96P 4 + 0.1P 5 + 0.32P 9 . 

Since /^ 3) = 1, /f } > 0, we have S h (X) > M 2 . Thus, fc(3) < 12.88 < 13. 

3-E. The one-sided kissing problem in four dimensions. Let H be a 

closed half-space of R". Suppose S is a unit sphere in H that touches the 
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supporting hyperplane of H. The one-sided kissing number B(n) is the maximal 
number of unit nonoverlapping spheres in H that can touch S. 

If nonoverlapping unit spheres kiss (touch) the unit sphere S in H C R™, 
then the set of kissing points is an arrangement on the closed hemisphere S+ 
of S such that the (Euclidean) distance between any two points is at least 1. 
So the one-sided kissing number problem can be stated in another way: How 
many points can be placed on the surface of S + so that the angular separation 
between any two points is at least 7r/3? In other words, B(n) is the maximal 
cardinality of a 7r/3-code on the hemisphere S + . 

Clearly, 5(2) = 4. It is not hard to prove that .8(3) = 9. Using extensions of 
Delsarte's method we proved that B(A) = 18 (see [25] for a proof and references). 
Recently have been obtained several new upper bounds for the one-sided kissing 
numbers 0[27l[3]. 

4 Spherical two-distance sets 

4-A. Two-distance sets. A set S in Euclidean space K" is called a two- 
distance set, if there are two distances c and d, and the distances between pairs 
of points of S are either c or d. If a two-distance set S lies in the unit sphere 
S™ -1 , then S is called spherical two- distance set. In other words, S is a set 
of unit vectors, there are two real numbers a and b, —1 < a, b < 1, and inner 
products of distinct vectors of S are either a or b. 

The ratios of distances of two-distance sets are quite restrictive. Namely, 
Larman, Rogers, and Seidel [T7] have proved the following fact: if the cardinality 
of a two-distance set S in R", with distances c and d, c < d, is greater than 
2n + 3, then the ratio c 2 /d 2 equals (k — l)/fc for an integer k with 

2<fc<i±^. 
~ ~ 2 

Einhorn and Schoenberg [13] proved that there are finitely many two-distance 
sets S in R ra with cardinality |5| > n + 2. Delsarte, Goethals, and Seidel [T2] 
proved that the largest cardinality of spherical two-distance sets in R™ (we de- 
note it by g{n)) is bounded by n(n + 3)/2, i.e., 

n(n + 3) 
g(n) < — - — . 

Moreover, they give examples of spherical two-distance sets with n(n + 3)/2 
points for n — 2, 6, 22. (Therefore, in these dimensions we have equality g(n) = 
n(n + 3)/2.) Blockhuis [7] showed that the cardinality of (Euclidean) two- 
distance sets in R™ does not exceed (n + l)(n + 2)/2. 

The standard unit vectors e\,... ,e n +i form an orthogonal basis of R" +1 . 
Denote by A„ the regular simplex with vertices 2ei, . . . , 2e n _|_i. Let A„ be the 
set of points ei + ej, 1 < i < j < n+1. Since A n lies in the hyperplane J2 x k =2, 



7 



we see that A n represents a spherical two-distance set in IR ra . On the other hand, 
A„ is the set of mid-points of the edges of A n . Thus, 

g(n) > |A„| - . 

For n < 7 the largest cardinality of Euclidean two-distance sets is g(n), where 
,g(2) = 5, g(3) = 6, g(4) = 10, g(5) = 16, and g(6) = 27 (see [H]). However, 
for n = 7, 8 Lisonek [21] discovered non-spherical maximal two-distance sets of 
the cardinality 29 and 45 respectively. 

4-B. Spherical two-distance sets with a + b > 0. In [29] using the polyno- 
mial method we proved the following fact: 

Theorem 4.1. Let S be a spherical two-distance set in R™ with inner products 
a and b. If a + b > 0, then 

\S\<^Lpl. 

Recently, Nozaki [3D] extended this theorem for spherical d-distance sets. 

Theorem 4.2. Let S be a spherical d-distance set in M. n with inner products 
ai, ■ ■ ■ , ad- Let 

n(t- flfc )=x:/ fe Gi" ) w- 

k=l k=0 

\s\< £ h k , 

k:f k >0 



Then 



where 

+ k-2\ /n + k-3 
k ) + { k-1 

4-C. Delsarte's method for spherical two-distance sets. Let S be a 

spherical two-distance set in M. n with inner products a and 6, where a > b. Let 
c = \J2 — 2a, d = v2 — 2b. Then c and d are the Euclidean distances of S. 
Let 

If k is defined by the equation: 6fc(a) = b, then (k—l)/k = c 2 /d 2 . Therefore, 
if |S| > 2n + 3, then k is an integer number and fee {2, . . . , K(n)} [17 . Here, 

Consider the case a + bk{a) < 0. Since bk(a) > —1, we have 



a £ Ik := 



2-fc 1 



2fe- 1 



Therefore, for fixed n, fc G {2, . . . , K (n)}, and a G we have spherical 
codes with two inner products a and (a) . The maximal cardinality of these 
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codes can be bounded by Delsarte's method (see details in [35] ). Since for 
6 < n < 40, n 7^ 22,23, Delsarte's bounds are not greater than n(n + l)/2, 
we have that g(n) = n(n + l)/2. For n — 23 we obtain g(23) < 277. But 
ff(23) > 276. This proves the following theorem: 

Theorem 4.3. // 6 < n< 22 or 23 < n < 40, then 

n(n + l) 
9{n) = 2 ' 

For n = 23 we have 

5 (23) = 276 or 277. 

The case n = 23 is very interesting. In this dimension the maximal number 
of equiangular lines (or equivalently, the maximal cardinality of a two-distance 
set with a + b = 0) is 276 18J. On the other hand, |A 23 | = 276. So in 23 
dimensions we have two very different two-distance sets with 276 points. 

Note that for n = 23 the Delsarte's bound is w 277.095. So this numer- 
ical bound is not far from 277. Perhaps stronger tools, such as semidefinite 
programming bounds, are needed here to prove that 9(23) = 276. 

Our numerical calculations show that the barrier n — 40 is in fact funda- 
mental: Delsarte's bounds are incapable of resolving the n > 40, k = 2 case. 
That means a new idea is required to deal with n > 40. 

5 Sylvester's theorem and SDP bounds for codes 

5-A. 2-point-homogeneous spaces. We say that a G-space M is a 2-point- 
homogeneous space if (i) M is a metric space with a distance p defined on it; 
(ii) M is a set on which a group G acts; (iii) p is strongly invariant under G, 
i.e. for x, x',y, 2/ G M with p(x,y) — p(x',y') there is an element g € G such 
that g(x) = x' and g(y) = y'. 

These assumptions are quite restrictive. In fact, if G is infinite and M is a 
compact space, then Wang 39J has proved that M is a sphere; real, complex 
or quaternionic projective space; or the Cayley projective plane. However, the 
finite 2-point-homogcneous spaces have note yet been completely classified (see 
for the most important examples and references [TU]). 

With any 2-point-homogeneous space M and an integer number k > are 
associated the zonal spherical function $>k(t) such that {<&fc(t)}fe=Q,i,2,... are or- 
thogonal on T := {t(x, y) : x,y £ M}, where r is the certain function in p 
( i.e r(x,y) = F(p(x,y))) defined by M. For all continuous compact M and 
for all currently known finite cases: <&&(£) is a polynomial of degree k. The 
normalization is given by the rule: $k{T~o) — lj where tq := t(x,x). Then 
$,,(*) = 1. 

Note that if M is a Hamming space F£ with r(x,y) — p(x, y)=Hamming 
distance, then To = 0. Here 3>fe(t) is the Krawtchouk polynomial Kk(t, n). Con- 
sider the case M = unit sphere §™ _1 C K" with r(x,y) — cos p(x,y) = (x,y), 
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where p(x, y) is the angular distance between x and y. Then the corresponding 
zonal spherical function $fc(t) is the Gegenbauer polynomial (t). 

5-B. The Bochner-Schoenberg theorem. The main property for zonal 
spherical functions is called "positive-definite degenerate kernels" or "p.d.k" 
[TU] , This property first was discovered by Bochner [5] (general spaces) and 
independently for spherical spaces by Schoenberg [3"4~] . 

Now we explain what the p.d.k. property means for finite subsets in M. 

Theorem 5.1 ([8, 34. 16J). Let M be a 2-point-homogeneous space. Then for 
any integer k > and for any finite C = {xi} C M the matrix ($fc(r(xi, Xj))) 
is positive semide finite. 

This theorem implies the fact that plays a crucial role for the linear pro- 
gramming bounds. For any positive semidefinite matrix the sum of its entries 
is nonnegative. Then 

Theorem 5.2 ( [TO US EE EE] ) • For any finite C = {x l } C M we have 

\G\ \C\ 

i=i j=i 

Since ^k{ T ( x i^i)) — < I ) fc(' r o) — lj this theorem implies 

i ^k(T(xi,Xj))>-l, k = 0,1,2,... . (5.1) 



5-C. The linear programming bounds. Let S be a fixed subset of T. We 
say that a finite subset C C M is an S-code if t(x, y) € S for all x,y £ C,x ^ y. 
The largest cardinality |C| of an S'-code will be denoted by A(M, S). 
The distance distribution {a t } of C is defined by 

at := (number of ordered pairs x,y £ C with t(x, y) =t). 

We obviously have 

«r = 1, «t > 0, t e T, ^ a t = (5.2) 

tST 

(5.1) and (5.2) make it possible to regard the problem of bounding A(M, S) 
as a linear programming problem: 

Primal problem (LPP): Choose a natural number s, a subset {ti, . . . , r s } of S, 
and real numbers a Tl , . . . , a Ts so as to 

maximize a T1 + . . . + a Ts 
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subject to 



a Ti >0,i = l,...,s, ^2 a Ti®k{n) > — 1, k = 0, 1, 



i=l 



This is a linear programming problem with perhaps infinitely many un- 
knowns at and constraints (1), (2). If C is an 5-code then its distance distribu- 
tion certainly satisfied the constraints (1), (2). So if the maximal value of the 
sum a T1 + ... + a Ts that can be attained is A*, then A(M, S) < 1 + A*. (The 



extra 1 arises because the term a r 



1 doesn't occur in this sum. 



Dual problem (LPD): Choose a natural number N and real numbers f±, . . . , /at 
so as to 

minimize /i + . . . + /jv 

subject to 

/ fc >o, k = i,...,N, 5^/ fc $fc(t)<-i, teS. 



k=l 



Thus, we have 



Theorem 5.3. If A* is the optimal solution to either of the primal or dual 
problems, then A(M, S) < 1 + A*. 

5-D. Sylvester's theorem. Sylvester's theorem (see [33J HZ]) gives an answer 
for the following question: Suppose we know for complex numbers t\,...,t n only 
its power sums Sfc, where Sk are real numbers. How to determine the number 
of real ti in an interval [a, 6]? 
Let 

/ s si . 
si s 2 . 



V 

/ si - as 
s 2 — as i 



- 1 Sm 

s 2 - as\ 
S3 - as 2 



S2m-2 / 



/ 6s - si 6s i - s 2 
6si - s 2 6s 2 - s 3 



\ 6s r , 



6S» 



'm+1 



S m GS m _i \ 
^m+l aS m 

S 2 m-1 — aS2m-2 / 

6s m _i s m 
6s m s m -]-i 



6s 2 -m.-5 



and 



fin 







H m (a,b) = H m (s ,si . . . ,s 2m -i, [a, 6]) := 



F+( ) 




S2m-1 / 
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Theorem 5.4. Consider real numbers t\, . . . ,t n in [a,b]. Then for any nat- 
ural number m the matrix H m (a,b) is positive semidefinite. Moreover, for 
m > rank(i?„) the converse holds, i.e. if H m (a,b) >z 0, then complex num- 
bers t\, . . . , t n with real Sk are real numbers in [a, b]. 

In fact, the constraint H m {a, b) y doesn't depend on n — sq. Indeed, let 

Sk ■= — — — — = — , H m (a,b) = -H m (a,b). 
n n n 

In other words, H m (a, b) can be obtained by substituting Sk for Sk in H m (a, b): 
H m (a,b) = H m (s~i, . . . ,s 2m ~i, [a,b]) := H m (l,s X: . . .,s 2m _i, [a, b]). 

Thus H m (a, b) h if and only if H m {a, b) t 0. 

5-E. Semidefinite programming. The standard form of the SDP problem is 
the following [37J EE|: 

Primal Problem: 

minimize c\X\ + . . . + cgxt 

subject to 

X h 0, where X = T lXl + ...+ T e x e - T a . 

Dual Problem: 

maximize (Tq,Y) 

subject to 

(T i ,Y)=c i , i = !,...£, 

y y o. 

Here To, T%, . . . , Tt, X, and Y are real Nx N symmetric matrices, (c%, . . . , q) 
is a cost vector, (xi, . . . , a^) is a variable vector, and by (A, B) we denote the 
inner product, i.e. (A, B) = Tr(AB) = J2 a ijbij- 

5-F. The SDP bounds. From here on we assume that is a polynomial 

of degree k, $fc(ro) = 1, and S = Tn[a,b] (the most interesting case for coding 
theory and sphere packings). 

In fact, the optimal solution A* of the LPP and LPD problems in 5-C 
depends only on the family of polynomials $ := {&k(t)}k=o,i,... and [a, b}. We 
denote 1 + A* by LP($, [a, 6]). 

Since 

we have 

fc 

$ k (t x ) + . . . + $ k (t e ) = J^PkdSd, s d = tf + ... + tf. 

d=0 
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Let C = {vi} be an S-code on M, and let Tjj = r(vi,Vj). Note that the 
number of ordered pairs (vi, Vj), i ^ j, equals n = \C\(\C\ — 1). Then (5.1) can 
be written in the form: 



k 



V + Pka + },PkdXd > 0, 



(5.3) 



where 



1 



X d = S d =—, S d = } T ir 



|C|-1 



From Theorem 5.4 we have for any m: 



H m (xi, . . .,x 2m -i, [a,b]) b 0. 



(5.4) 



Now we introduce the simplest SDP bound. 

SDPo Problem: Choose a natural number m and real numbers y, xi, . . . , £2m-i 
so as to 



Note that in (5.3): |C| = (1 + y)/y. Thus 
Theorem 5.5. If y* is the optimal solution of the SDPq problem, then 

A(M,S) < SDP (<S>,la,b},m) := l±f-. 

y 

Since we just substituted H m y for t € S in the LPP problem, we can 
expect that S , L>P (*, [a, 6]) = LP($, [a, b}) (see details in [27]). 

In fact, for a continuous M the LPP and LPD problems are not finite linear 
programming problems. These problems can be solved only via discretization. 
For instance, Odlyzko and Sloane [31] (= [TOj Chapter 13]) applied LPD for 
upper bounds on kissing numbers, where they replaced S by 1001 equidistant 
points in S. For the LPP problem it is not clear how to do a discretization 
of {a T }. On the other hand, for a given m the SDPq is a finite primal SDP 
problem. As a by-product of solutions of this problem we have bounds on \C\ 
and power sums s k (see [23 Section 5]). 

In [571 Section 6] is shown that some recent extensions of Delsarte's method 
can be reformulated as SDP problems (SDP A). Section 7 extends the SDP A 
bounds to subsets of a 2-point-homogcncous space and shows that some upper 
bounds for codes can be improved. In particular we obtain new bounds of 
one-sided kissing numbers. 



minimize y 



subject to 



k 




1, ... ,2m - 1, 



H m (xi, . . .,x 2m -i, [a,b]) b 0. 
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6 Multivariate positive definite functions 



6-A. Schrijver's method. Recently, Schrijver [35] using semidefinite program- 
ming (SDP) improved some upper bounds on binary codes. Even more recently, 
Schrijver's method has been adapted for non-binary codes (Gijswijt, Schrijver, 
and Tanaka [TS]), and for spherical codes (Bachoc and Vallentin [51 [3113]). In 
fact, this method using the stabilizer subgroup of the isometry group derives 
new positive semidefinite constraints which are stronger than linear inequali- 
ties in the Delsarte linear programming method. We consider and extend this 
method for spherical codes in [38J . Note that this approach is different from the 
method considered in Section 5. 

6-B. Multivariate positive definite functions on spheres. Let us consider 
the following problem: For given points Q — {qi, . . . , q m } in M to describe the 
class of continuous functions F(t, u, v) in 2m + 1 variables with t€R, u, v € 
R m , F{t, u, v) = F{t,w,\i) such that for arbitrary points p\, . . . ,p r in M the 
matrix (F(tij,Ui,Uj)) hO, where Uj = r(pi,Pj), Uj = (r(pi, q±), . . . , r(pi, q m )). 

Denote this class by PD(M, Q). If Q = 0, then PD(M, Q) is the class of p.d. 
functions in M. In this case an answer is given by the Bochner - Schoenberg 
theorem. 

Let < m < n - 2, t € R, u, v e R m for m > 0, and u = v = for 
m = 0. Then the following polynomial in 2m + 1 variables of degree k in t is 
well defined: 

4 n ' m) (i,u,v) := (1- |u| 2 ) fc / 2 (l- |v| 2 ) fe / 2 Gi"- m) ( *-< u ' v > ^ 



Vv/(i-H 2 )(i-ivi 2 )y 

In [55] we proved the following theorem: 

Theorem 6.1. Let < m < n - 2. Let Q = {qi,...,q m } C § n_1 with 
rank(Q) = m. Let e±, . . . , e rn be an orthonormal basis of the linear space with 
the basis q\,..., q m , and let Lq denotes the linear transformation of coordinates. 
Then F G PD(§"-\ Q) if and only if 

OO 

F(t, u,v) = J2 ^ v ) G i"'" l) (*' L Q (u),L Q {v)), 

k=0 

where /fe(u,v) y (i.e. /fe(u,v) = /ii(u)/ii(v) + . . . + hg(u)hi(v)) for all k > 0. 

Note that for the case m = that is Schoenberg's theorem [Mj, and for 
m = 1 it is the Bachoc- Vallentin theorem [3]. 

6-C. Upper bounds for spherical codes. In this subsection we set up upper 
bounds for spherical codes which are based on multivariate p.d. functions. These 
bounds extend the famous Delsarte's bound. Note that for the case m = 1 this 
bound is the Bachoc - Vallentin bound [3]- 
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Definition 6.1. Consider a vector J = (ji, . . . ,jd)- Split the set of numbers 
{ji, . . . ,jd} into maximal subsets I\, . . . , with equal elements. That means, 
if I r — {j ri , ■ ■ ■ ,jr s }, then j ri = . .. = j rs — a r and all other je ^ a r . Without 
loss of generality it can be assumed that i\ = \Ii\ > . . . > = \I k \ > 0. (Note 
that we have i\ + . . . + = d.) Denote by ip(J) the vector u = (ii, . . . ,ife). 
Let 

W d := {u> = (ii, . . . , i k ) : i\ + ■ ■ ■ + ih = d, i\ > . . . > ik > 0, h, . . . , i k £ Z}. 
Let w £ Wd- Denote 

UN) := #{J = (ii, {1 N} d : ^(J) = w }, 

UN) 



qUN) :-- 



N 



It is not hard to see that if u> £ Wd, then qu(N) is a polynomial of degree 
d—1, and 

]T qUN) = N d -\ 

Definition 6.2. For any vector x = with 1 < i < j < d denote by A(x) 

a symmetric d x d matrix (ajj) with all 
Let < 9 < 7T and 

X(0) := {x = {xij} : Xij £ [— l,cos6>] or Xij = 1, 1 < i < j < d}. 

Now for any x = {:%} £ X(9) we define a vector J(x) = (ji, ■ ■ ■ ,jd) such 
that jk = k if there are no i < k with x ik — 1, otherwise ju — i, where i is the 
minimum index with Xik = 1. 

Let u) £ Wd- Denote 

D u (0) := {x £ X{6) : ^(J(x)) = 10 and A(x) h 0}. 
Let /(x) be a real function in x, and let 

B u (0,f):= sup /(x). 

Note that the assumption A(x) y implies existence of unit vectors pi, . . . ,p d 
such that ^4(x) is the Gram matrix of these vectors, i.e. Sy = (pi,pj). More- 
over, if = 1, then pi = pj. In particular, D^ d ){9) = {(1, . . . , 1)} and therefore 
B (d )(0,/) = /(l,...,l). 

Definition 6.3. Let x = {s t3 }, where 1 < i < j < m + 2 < n, and let A(x) y 0. 
Then there exist P = {pi, ■ • • , Pm+2} C § n_1 such that x^ — (pi,Pj)- Let 
F(x) be a continuous function in x with F(Siki) = Fix) for all x^i that can be 
obtained by interchanging two points pu and pi in P. We say that F(x) £ PDJ^ 
if for all x with A(x) y we have F(xi 2 , Ui, u 2 ) £ PD(S" _1 , Q(xj), where 
Ui = (x i3 , ■ ■ .,Xi, m +2), Q(x) = {p 3 , ■ ■ ■ ,Pm+2>, and F(x 12 ,u 1 ,u 2 ) = F(x). 
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For the classical case to = Schoenberg's theorem says that / € PDq ^ an d 
only if 

with all fk > 0. It is not hard using Theorem 6.1 to describe class of functions 
in PD^ for all m < n - 2. 

Let C be an iV-element subset of the unit sphere S n_1 C R n . It is called 
an (n, N, 9) spherical code if every pair of distinct points (c, c') of C have inner 
product (c, c') at most cos0. 

In [SB] we proved the following theorem. 

Theorem 6.2. Let f > 0, < to < n - 2, and F(x) = /(x) - / € PD™ . 
XTien an (n, AT, 9) spherical code satisfies 

foN m+1 < J2 B u (0,f)q u (N). 

It is easy to see for m = that 0( 2 )(A r ) = 1, 9(1,1) {N) = N — 1, and 
B(2){9,f) — /(l). Therefore, from Theorem 6.2 we have 

/ iV</(l) + S (1 ,i)(^/)(iV-l). 

Suppose B (1A) (9J) < 0, i.e. /(*) < for all t e [-l,cos0]. Thus for (n,N,6) 
spherical code we obtain Delsarte's bound 

fo 

The Bachoc-Vallentin bound [H Theorem 4.1] is the bound in Theorem 6.2 
for to = 1 and B {1 1 1} (0, /) < 0. Indeed, let B (2 ^(0, /) < B. Since q {3) (N) = 
1, q (2A) (N) = 3(N- 1), and B {3) (9,f) = /(l, 1, 1), we have 

JoiV 2 </(l,l,l) + 3(7V-l)S. 

Let us consider Theorem 6.2 also for the case m = 2 with i?(i i i n(0, /) < 0. 
Let B( 3 ,i)(0,/) < Bi, B( 2l 2)(0,/) < B 2) and B (2 ,i,i)(0, /) < B 3 - Then 

/oiV 3 < /(1, 1, 1, 1, 1, 1) + 4(JV - + 3(7V - 1)B 2 + 6(JV - l)(N - 2)B 3 . 

Let /(x) be a polynomial of degree d. Then the assumptions in Theorem 
6.2 can be written as positive semidefinite constraints for the coefficients of F 
(see for details [3JGI1IUE3GIS]). Actually, the bound given by Theorem 6.1 can 
be obtained as a solution of an SDP (semidefinite programming) optimization 
problem. In [2l [3] using numerical solutions of the SDP problem for the case 
m = 1 has obtained new upper bounds for the kissing numbers and for the 
one-sided kissing numbers in several dimensions n < 10. 

However, the dimension of the corresponding SDP problem is growth so fast 
whenever d and m are increasing that this problem can be treated numerically 
only for relatively small d and small m. It is an interesting problem to find 
(explicitly) suitable polynomials F for Theorem 6.2 and using it to obtain new 
bounds for spherical codes. 
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